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XX. Coordinate Transformation 


In order to facilitate the coincidence calculations, the coor¬ 
dinates of the ship and the satellite are transformed to the ECO 
system, in which the equatorial plane is the plane of the satellite's 
orbit. The transformation from earth-centered fixed (ECP) to ECO 
coordinates is shown in Figure 1. 

These are the transformation matrices for each step: 



where 

(2) “t = “o + e rot* (t ■ V 

a Q is the Greenwich hour angle (GHA) of y , the first point of 
Aries, at t 

o 

e rot eartia ' s rotation rate 

Noto: a t , measured in radians, is the Greenwich sidereal time at 


Universal Time t. 







Step (1) 


Step (2) 
ing node 


step (3) 


ECP to ECI. 


z,z 

A. 


Rotate the equatorial plane 
by the Greenwich hour angle 
at time t (U.T.),* a*.' 


■ > y 


Rotate ECI coordinates by ft, the longitude of the ascend- 



Tilt the equatorial plane by the inclination of the orbit. 


Z* 



Figure 1 
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cosft 


sinft 0 


sinft cosft 0 


(4) A « [ 0 cos i sin i 

\ 0 -sin i cos # 


The complete transformation is: 


cos (ft - O 


(J5) ^ - 2L A A = ( (-cos i) sin (ft - ct^) (cos i) cos (ft - a^) sin i 

\(sin i) sin (ft - a fc ) (-sin i) cos (ft - a fc ) cos i 


The ECO coordinates of the ship are: 


(6) i Y or k = cos i [-X sin (ft - a fc ) + Y cos (ft - a fc ) ] + 2sin i 
z orb “ s i n (^ “ ” Y cos (ft - a^.)] + Zcos i 


where X, Y, Z are the ship's ECF coordinates. 
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XIX. Calculation of Time Intervals 

The ship can be observed when it is in a band about the equator 
in the ECO system. The width of the band U lim ) Is determined by the 
scan pattern of the instrument. The first step in obtaining the times 
the ship enters and leaves the band is to set 2 orb * + in (6) and 

solve quadratic equations for cc?r ct fc and sin a fc (see Appendix 1) . We 
need both sin a fc and cos x { , so that we can put in the proper qua¬ 
drant. We then use (2) to calculate the U.T. times the ship is at the 
borders of the band in an Interval of approximately 48 hours centered 
about the time of the ship's observation. Depending on the ship's 
latitude, there will be two, three, or four time intervals to be 
searched for coincidences. 
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IV, Existence and Times of Coincidence 

The work of OTIME is nearly complete after the time intervals 
have been calculated. It calls LONGIT, which checks for the exist¬ 
ence of a coincidence in each interval using the criterion; 

(lon sat (t 2> " lon ahlp <V> * 0 (7 >' 

If the time interval is longer than the period of the satellite, 

LGNGXT is restricted to searching the full period closest to the 
time of ship observation. Furthermore, if the interval is longer 
than half a period, it is split into three parts in order to avoid 
a breakdown of (7). 

GETTIM uses a modified regula falsi technique to find the time 
of coincidence. The first step in the procedure uses the normal 
regula falsi technique, The number of iterations to be done is 
passed to the coincidence program from the main program. If conver¬ 
gence has not been reached after the iterations have been done, the 
second step is to determine in which half of the interval the root 
lies and repeat the first step using that half-interval. The max¬ 
imum number of times the interval can be halved is also passed from 
the main program. 

If convergence is reached, an empirical correction factor, which 
depends on the month being processed, is applied to the time of coin¬ 
cidence. The final time of coincidence is 


where t Q = originally calculated time of coincidence (in hours) 
Wrl' t corr2 = correct i°n factors. 


( 8 ) 


The correction factors are obtained from plotting <5t, the error in 
the calculated SMMR observation time, against At, the difference 
between the time of the ship observation and the SMMR observation time, 
Points which were observed by SMMR are read from the SMMR TETs for 
three or four days spread over the month being processed, These 
points are fed to OTIME along with several ship observation times. 
These points fit a straight line whose slope, in minutes per hour, is 
^corrl and whose fit-intercept is the negative of t corr2 * 



The following shows how to us© OTIME and describes the output. 

Unless otherwise noted, all variables listed below are HEAL*8. 

A. Input to OTIME 

ALON, ALAT The position of the ship. Units are degrees. 

The range of ALON is (0, 360), and it is 
measured positive to the west . The range of 
ALAT is (-90, +90). 

TLIM In GETTIM, the interval of coincidence is 

reduced until its length is leas than the 
value of TLIM. A value of ,01 minute is 
recommended. 

NITER(I*4) Also used in GETTIM, NITER is the number of 
iterations to be done before the time inter¬ 
val is halved. A value of 3 .is recommended. 

NHALF(1*4) The maximum number of iterations is the pro¬ 
duct of NHALF and NITER. A value of 17 is 
recommended. 

ITSOBS(1*4) The hour (U.T.) of the day at which the ship 
made the observation, ITSOBS must be in the 
interval (0,24), 

IDAYYR(1*4) The day of the year on which the ship made 

the observation. Midnight of January 1, 1978 
is day 0.0. 

LYEAR(1*4) The year in which the observation was made, 
LYEAR must have 4 digits, i.e., 1978 not 78. 



TC0RR1 


Time correction factora used in GETTXM. 

Units are minutes/hour, TCORR1 is read in 
by the namelist THIS from unit 9. 

TC0RR2 Time correction factor used in GETTXM. Units 

are minutes, TC0RR2 is read by the namelist 
TIME from unit 9, 

When OTXME is called for the first time, it reads in TCORR1 and 
TC0RR2 from unit 9, and from unit 11 it rsads in a set of orbital 
elements and other information which enables it to calculate the 
ECX coordinates of the satellite at any time during the month which 
is being processed. The following cards must be included in the JCL 
for the GO steps 

//GO.FTQ9F0Q1 DD * 
namelist input' 

and 

//GO-.FT11F001 DD DSN-2MMHT. ELTS, aaann, DXSP*SHR 

where aaa » first three letters of the month being processed 
nn * last two digits of the year 

B, Output from OTXME 

TOBS: 4-element time vector containing times of coinci¬ 
dences measured relative to the time of observation 
(l.e., t obs 85 0). Units are hours (U.T.) . 

IFOBS: 4-element flag vector. If IFOBS(n) » 0, no coin¬ 
cidence occurred, but if IFOBS(n) =* 1, a coincidence 
occurred, and TOBS(n) contains the time. 


IERFs 


Error flag (see Sec, C). 




C. Error Conditions 


There is on® «rror flag, IERE, 
Valu® of t erf 

a. Error conditions sat in OTIHE 
0 
1 


which can tak® several valuess 
Condition 

Normal 

The elements passed to GKEP 
violated one of the follow¬ 
ing conditions? 

(i) The eccentricity was 
not in the range (0, 1). 

(ii) The semi-major axis 
was not in the range 
(0, 1 x 10 35 ). 

The year in which the 
observation was made does 
not match the year for 
which the orbital elements 
were read in. 


3 


4 

5 

6 

7 


Orbital elements were read 
in for a range of days 
which does not include the 
day the observation was 
made on. 

The hour of the day is out 
of range. 

The latitude of the ship 
is out of range. 

The longitude of the ship 
is out of range. 

The absolute value of a 
computed time is greater 
than 35 hours. 
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8 The times that the ship 
entered the band is greater 
than the time it left the 
band. 

9 The ship remained in the 
band for more than one 
sidereal day, 

15 One or both of LDAYQ and 
&DAYF, parameters which 
point to the first or last 
days for which orbital 
elements have been calcu¬ 
lated, are out of range, 
or LDAYO is greater than 
LDAYF. 

16 Orbital elements do not 
exist for the day given by 
IDAYYR 


b. Error conditions set in LONGIT. 


10 

11 

12 


Same as 7. 
Same as 8. 
Same as 9. 


c. Error conditions set in GETTIM. 

13 The differences in long¬ 
itude between the ship and 
the satellite have the same 
sign at both endpoints of 
the time interval. 

14 The time of coincidence 
(U.T.) has an absolute 
value greater than 30 
hours. 



NHALF*NITER Convergence was not reached, 

If an error Is encountered, an error message is printed, control 
is returned to the main program, and no ruther processing is done. 


D. Common Block 

The common block CONST appears in all subroutines except GKEP, 
The following table lists the variables contained in CONST. They 
are all computed in OTIME; 


X/Y, z 
ALONO 

SREVH 

COSI, SINI 
Cl, C2 

T 


ECP coordinates of ship. 

Longitude of satellite in ECO system at time 
of ship observation. 

Revolution rate of satellite about the earth 
in radians per hour (U.T.). 

Cosine and sine of orbital inclination. 
Constants used in calculation of X ^ and 
Y orb ^ see f and a t ( See App * i) ' 

Eight element vector containing endpoints of 
time intervals to be searched for coincidences. 


TP I 2 tt . 

EROTH Rotation rate of earth, e rofc , in radians per 

hour (U.T.). 


GHAO 

TPIE 

SINOM, COSOM 

ERSINI 

PERIOD 


GHA at time of ship observation. 


Sine and cosine of the longitude of the 
ascending node. 


_e rot ' sin (i) where sin(i) = SINI. 
Period of satellite orbit as computed by 
Kepler's third law. 


NTIMES The number of time intervals to be searched 

for coincidences is NTIMES divided by 2. 
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E. Calling Sequence 

CALL OTIME(ALON, ALAT, TOBS, TLIM, IFOBS, IERF, NHALF, NITER, 
ITSOBS, IDAYYR, LYEAR), 


Appendix 1. Calculation of 

(Al-1) ^ ±z ii ra / s ^ n *•) “ [ x sinQ <- y cosft] cos « fc 

- [x cosfl + y sinQ] sin + z cot i 


x sinft - y cosft 
x cosQ + y sinfi 


(Al-3) 


sin a. =Ji~^ 


cos a, , then (Al-2) becomes 


B 2 (1 - cos 2 a, ) = A 2 cos 2 a, + 2AC cosa, + C 2 






Appendix 2. Block Diagram 
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APPENDIX 3 
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:r.sc, 7/7! 




ATX Xb;*D/0/,GB/.3££6032D6/ 

aTA ZllB/?lOO.B3/,3in:»Y/23.93tja695C3/ 

ATX PI/3. 1*11592£53565292D0/,R i/6.3G71X25P3/ 



b-2(H,1)=Sr.BI-HXJ09 AXIS (KILOHKT.lnSJ 








Uva/SKfiatH) hoiiou KiLa=(9'*) 



It-1 (3) = ’.ITS’ (IBHY.3) 






C050H=CC05 (X*CI («1) 1 

GUXO=EBOC ((5L152 (ICAI.9) »6.3303B736DO*DTSOES) ,TPIJ 
TH !T*0=X’CI 








l 3=Z/MJUJ ( t iITS2 (1C AY , 3) ) 




CjiLCUUTT I1BJ3 CH INriFYAL 0,21) 














OZIH=DZCRET (K*1) 










1/0 IF |T{K).GT.T {»♦!)) T (K) =1 (K) -TPIE 



0) HSUJHH 




















> V 





UEEOUHKS lCKGII (IfCES.IiBF) 



X’i= (TOC1-HN) -(EEItIOt/2.C0) 




IP ( (CABS (ICB'f-TJH) .Gl.EZHIOD) .IKD. (DAaS(TIK) .GT.DUBS ITOOT))) 



J Of 


IFC-IjI |U + 1}/2} = 1 
1 (U)=llnc (2»IFI2-1) 




XF C (CAES (TiK) ) .Gl. 35.DC) .OB. (DABS (X (K*1) ) . GT.35.P0) ) IlsRF=10 
xF(IXHF.M.O) MBITS (6,1) 

XF(1BBF.B>.0) 8Z1UFN 














=02, LI*eCHT=56 , SIZli= 3000!? 






s*na= BiiN,cFX=o?,i.m%csx= 56 ,siri=ooooK, 

SOUHCa^BECDIC^KCSS^ii'.HOmCK.LOID.JUP.Jimnxr.ID.XKrtP 





) 0 













SOUBC'k, 1EODIC,KCII3 f,]IODuCK ,LOUD,HHP, NOi’I)IT,ID,Ifii;F 
LOODLfe PH?CJ31CH FUkCTJCK TORB (GEAJ 



) 



roaB / 



) # 
















OOOOOOOaQOOOOOOOeOOOOOoCOOOOOQOQOQQOOOOOOQOOOOOOOOOOOO 





iF ( (»L*3E1.n.PI) - UIS^Ta.t^.TFI) ) 4I.H*5T1 =M.HST 1 -3FI 













I 

1 

1 













LI «3CSI=56 ,SIZl*=OOO0 K 








El rt 79.205/15.sa.no 




















• CJD # HSE # MOWIT,Ii>,X*l>F 





